Generalized problem of two and four Newtonian 



centers 

Alexey V. Borisov, Ivan S. Mamaev 

Institute of Computer Science, Udmurt State University 

1, Universitetskaya str., 426034 Izhevsk, Russia 

Phone/Fax: +7-3412-500295, E-mail: borisov@rcd.ru, mamaev@rcd.r 

February 8, 2008 

Abstract 

We consider integrable spherical analogue of the Darboux potential, which 
appear in the problem (and its generalizations) of the planar motion of a 
particle in the field of two and four fixed Newtonian centers. The obtained 
results can be useful when constructing a theory of motion of satellites in the 
field of an oblate spheroid in constant curvature spaces. 
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1 Introduction. 

The two-center problem is well known in classical celestial mechanics: two fixed 
centers, with masses m l and m 2 , attract some "massless" particle, moving in their 
field according to Newton's law. The integrability of this problem was proved by 
Euler, by means of the separation of variables j!5j . 

A qualitative analysis of the plane two-center problem was offered by C. Char- 
lier [7] (see also |H]); a qualitative analysis of the spatial two-center problem can be 
found in the paper by V. M. Alekseev Note also that it was Lagrange who ob- 
served that the two-center problem remains integrable upon addition of the potential 
of an elastic spring attached to the midpoint of the rectilinear segment connecting 
both centers. Lagrange also studied the limiting case of the problem, where one of 
the two-centers and its mass tend to infinity. In the limit we have the problem of 
a point subject to the superposition of the field of a Newtonian center (the Kepler 
problem) and a homogeneous field. The corresponding separation of variables and 
a qualitative study of this problem was done by M. Born in his book on atomic 
mechanics jU], in connection with the splitting of a hydrogen atom's spectral lines 
observed when it is put into an electric field (the Stark effect). 

A more general two-dimensional integrable system, which incorporates the two- 
center problem, was found by G. Darboux (1901) [10J, by means of the method of the 
separation of variables. In this paper Darboux also obtained the existence conditions 
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for an additional quadratic integral in case of a natural system on a plane. Later, 
these conditions were described by Whittaker fZT\ . 

Consider a particle of unit mass, moving on M 2 = {x, y} in the potential field 

A A' B B' B l B[ „ 2 

v = - + — +-+—+— + ^ + Cp\ (1) 

x y r r r 1 r\ 
where A, A', B, B', B 1 , B[, C = const. Here, r and r' are the real distances 
between m and the two identical real centers positioned at the points (— c, 0), 
(c, 0) on the abscissa axis, r = ^JJx — c) 2 + y 2 , r' = + c) 2 + y 2 \ p is the 

distance between m and the origin 0\ r 1 and r[ are the "complex distances" be- 
tween m and the two imaginary centers (0, di) and (0, — di), r x = x 2 + (y — id) 2 , 



r[ = a/x 2 + (y + id) 2 (see Fig. 1). 

For the potential (0) to be real, it is necessary that B[ and B\ are complex 
conjugate: B\ — B\. As is shown in ^Oj, if d = c, the system (0) allows the 
separation of variables in elliptic coordinates 

x = c chv cos u, y = c shv sin u 

and has an additional first integral, which is quadratic in the momenta. 

Let us examine some special cases of the potential The case of B 1 = B[ = 
was studied by G. Liouville (a more special case of A = A' = B l = B[ = was, as 
we already mentioned, described by Lagrange). 

It was shown in the paper pQ that the problem of a particle moving in the field 
of two complex conjugate centers, i. e., when A = A' = B = B' = C = in (0), is 



integrable in the three-dimensional space and makes a good approximation to the 
problem of a satellite's motion in the field of an oblate spheroid (eg., the motion of 
an artificial satellite of the Earth). 

In the paper by I. S.Kozlov ^7], the problem of the plane motion of a parti- 
cle in the field of four fixed centers (two real and two complex) was integrated in 
terms of quadratures and further studied. Besides, in [T7j, several interpretations 
of this problem were offered, with reference to actual problems of applied celestial 
mechanics. 

2 The Kepler problem. The two-center problem 
on a sphere and a pseudosphere. Historical 
notes 

Systematic generalization of various problems from classical and celestial mechanics 
to constant curvature spaces (a three-dimensional sphere S 3 , as well as a pseudo- 
sphere L 3 , or Lobachevsky space) was done by W. Killing in his extensive, but, 
unfortunately, almost forgotten paper [T^j . 

Note also that beside Killing, in the 19th century, non-Euclidean mechanics in 
constant curvature spaces was studied by R. Lipschitz, F. Schering and H. Liebmann. 
It is interesting that though a whole chapter from Liebmann's textbook on non- 
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Euclidean geometry (20] concerned the generalization of Newton's law of attraction, 
study of the Kepler problem, and reformulation of Kepler's laws for the cases of a 
sphere and a pseudosphere, similar results were independently and almost simultane- 
ously rediscovered in the 20th century by several authors [T^llSl ITBllTHll^l ITU 1231 IU] . 
The classical paper by E. Schrodinger [22] should also be mentioned, where he stud- 
ied a quantum analog of the Kepler problem in a curved space, implicitly assuming 
that the corresponding classical problem was integrable. An analogue of Newton's 
law of attraction for L 3 was known to J. Bolyai, and N. I. Lobachevsky, and for S 3 - 
to P. Serret. 

W. Killing in ^1 studied, among other things, the problems of n-dimensional 
dynamics in constant curvature spaces, including the dynamics of an n-dimensional 
rigid body. An up-to-date analysis can be found in ^2] (see also [5]). 

The generalization of the two-center problem to constant curvature spaces is also 
due to W. Killing, who integrated this problem, using the method of the separation 
of variables. It was independently solved in [T^l, where a more general problem 
was studied, similarly to what Lagrange did by introducing an elastic interaction 
potential in the plane two-center problem. In the papers j2El I2H] , a bifurcational 
analysis of the two-center problem on a sphere and on the Lobachevsky plane was 
offered. In [5] we examined the spatial two-center problem from the standpoint of 
its reduction and integrability; we also studied other integrable and non-integrable 
problems of celestial mechanics in curved spaces (including the restricted two- and 
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three-body problems, behaviour of libration points, dynamics of rigid bodies). 

In this paper, we will offer an explicit algebraic expression for the first integral 
of the generalized two-center problem from ff3 [I^j, and show a new analogue of 
the problem of four Newtonian centers and n Hookian centers. In this paper, we 
study only the case of a two-dimensional sphere S 2 , though all the reasoning can 
easily be extended to a pseudosphere L 2 . Certain (not all) results are generalized 
to a three-dimensional sphere S 3 (or a pseudosphere L 3 ). 

3 Generalization of the two-center problem to S 2 . 
Additional quadratic integral 

Assume that a unit sphere S 2 is given in the three-dimensional space M. 3 = {q l , q 2 , q 3 } 
by |q| 2 = q\ + ql + q% = 1 and denote by q = (q 1 , q 2 , q 3 ), p = (p la p 2 , p 3 ) the 
redundant coordinates and momenta, respectively. Now if we introduce the an- 
gular momentum vector M = p x q and put 7 = q, it is easy to show El El 
that the equations of motion in an arbitrary potential V = V(q) = V(-y) can be 
presented as a Hamiltonian system with the Poisson bracket defined by the alge- 
bra e(3) = so(3) ® s M 3 : 

{M„ Mj} = e ijk M k , {M { , Mj} = e ijk M k , { 7i , 7j } = (2) 
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and the Hamiltonian 

H=^M,M) + V(j). (3) 
From (J2J), (J3J) we obtain the equations 

dV 

M = 7 x — — 7 = 7 x M, 
cry 

which coincide with the equations of motion of a spherical top in the potential V("y) [HI 
0. 

The bracket (j2J) is degenerate and has two Casimir functions: F 1 = (M, 7), 
F 2 = (7, 7) = 1. For the problem of a point moving on a sphere, it is necessary 
that F x = (M, 7) = (p x 7, 7) = 0. 

It is well known that the analogues of the Newtonian and Hookian potentials 
on S 2 are, respectively, U l = ficotO and U 2 = ctan 2 #, fi, c = const, where 9 is 
measured from a certain fixed pole on the sphere [T6| I19j. 

Consider the potential 

V = cot 9 X — ji 2 cot 9 2 , (4) 

where the intensities of the Newtonian centers, while 9 i is the angle 

between the radius- vector of a particle and the radius- vector of the z-th center. Place 
the Newtonian centers at the points r 1 = (0, a, f3),r 2 — (0, —a, /3), a 2 +/3 2 = 1, and 
add to (HJ), for the sake of generality, the potentials of three Hookian centers located 
at the mutually perpendicular axes \ cj^f (Cj = const). We should also introduce 
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the quadratic potential C(a 2 7 2 — /3 2 7 2 ) C ^ 0, which is a special case of Neumann's 
potential. On the level (M, 7) = we find two commuting functions {H, F} = 0, 
quadratic in M pi2"T|: 



P>3 ~ a l2 , 1 72+73 



\/7i 2 + /9i72 2 + « 2 7l + 2a/3727 3 2 ' 7? 
1 7i 2 + 7| , 1 7i + 7 2 2 , r( 2 2 ~2 2X (5) 

+ o C 2 To +o C 3 -5 + 72 -P 7 3 J, W 



F = a 2 M 2 - (3 2 M 2 + 2a(3{V 1 - V 2 )- 



— \P 7 2 - « 7 3 ) 2P 7i + — a 7i + 2Ca P 7u 

7i 72 7s 



where /i 2 , a, /3, c 1; c 2 , c 3 , C = const, and the functions V^, V 2 are: 

/i 1 (/?7 2 + «7 3 ) 



*2 



V / 7i 2 + /5 2 72 2 + « 2 73 2 - 2«/3727 3 : 
/i 2 (p7 2 - «7 3 ) 



(6) 



+ /3 2 7 2 2 + a 2 7 3 2 + 2ap7 2 7 3 ' 
The function if is the Hamiltonian, and F is an additional quadratic integral. As it 

is noted in j^j, the integrability of the corresponding three-dimensional (S 3 ) problem 

closely depends on whether a Hookian center (with potential c/7 2 ) can be added (at 

some point of the arc joining the Newtonian centers) to the two-center problem 

without violating the problem's integrability. Indeed, the term c/7 2 , c = const 

appears in the three-dimensional CclSG clS 8b result of the Routh reduction procedure, 

which uses the cyclic integral. This integral is due to the equations' invariance under 

rotations (group 50(2)), in the plane perpendicular to the plane of the two centers. 
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The system (J5J) is of the Liouville type and can be integrated in sphero-conical 
coordinates u 1 , u 2 , (0 < u Y < a, < u 2 < (3) given by 

7 2 = a/ (a 2 — u 1 )(a 2 + u 2 )/a, (?) 

7 3 = ^((3 2 + Ul )((3 2 -u 2 )/(3. 

Note, however, that obtaining the integrals © in the algebraic form is quite a 

non-trivial problem, as its solution implies dealing with an inverse sphero-conical 
transformation. 

As A. Albouy informed us, the two-center problem on S 2 (or L 2 ) can be trans- 
formed to the traditional Euler problem of two centers by means of the central 
(gnomonic) projection and a suitable transformation of time. However, we cannot 
yet prove this statement. 

4 The problem of four Newtonian centers on S* 2 

Consider the potential on a sphere: 

Vjm = ^l COt^+^COt^ = 

= ^ M7i±!^3 + ^ M7i ~ ij^h ^ 

yj (fi 2 - v 2 ) 2 - (>7i + ^7 3 ) 2 V (M 2 - ^ 2 ) 2 - 07i - ^7 3 ) 2 ' 

where /i 2 — z/ 2 = 1, £ 1; £ 2 = const. 

This potential corresponds to the two-center problem on sphere. The intensities 
of the centers are "complex" and the centers themselves are equidistant (complex 
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conjugate) from the pole (Fig. 2). For the potential to be real, it is necessary 
that £ x = £ 2 . As in the Euclidean case, the potential (jHJ) can be regarded as a 
certain approximation to the problem of a particle moving in the field of an oblate 
spheroid in a curved space. 

The system with the potential (jHJ) is also separable in the sphero-conical coordi- 
nates ((7j), provided that 



(3 af3 
1 — cr I — a 2 



In the coordinates (J7|), the potential V + V Im is also separable. This potential 
(for Cj = 0) corresponds to the problem of four fixed centers, two imaginary and 
two real, which belong to two mutually perpendicular planes through the pole (see 
Fig. 2). Here, as in the planar case, when the distance between the real centers is 
fixed, the distance between the complex centers is also not arbitrary: it is uniquely 
defined by ©. 

It is easy to show that the potentials 

can be added (without violating integrability) to the four centers problem and it 
results in a more general system, separable in the coordinates ((7j). The poten- 
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tials V, V lm , V G , V n written in terms of the variables (J7|) look like: 

Oi + l* 2 ) a/(« 2 - u^iP 2 + u x ) + - fi 2 ) a/(« 2 + u 2 )((3 2 - u 2 ) 



V 



Ui + u 2 

71 % + M 2 

1 2 (a 2 -« 1 ) _1 -(a 2 + u 2 ) _1 

2 ^ I ' 

72 M l + m 2 

— o = ap ■ , 1/^ = — 



73 ' M l + ^2 " «1 + «2 

One can easily show that in the limit R — > oo (the case of Euclidean plane) the total 
potential V + V| m + V G + V N becomes the Darboux potential ((H). Note that this po- 
tential, or even V + V^, can no longer be generalized to the corresponding integrable 
potential of the three-dimensional problem (S 3 ), because there is no cyclic integral, 
though, taken individually, the potentials V and V\ m allow such a generalization. 



5 The problem of n Hookian centers on a sphere 

Let us present one more integrable modification of the problem of a mass point 
moving in the field of the Hookian potentials 0^/(7, rj 2 , c i = const, where the 
Hookian centers of attraction r t ,i = 1, 2, . . . , n do not belong to mutually orthogonal 
axes, but are placed arbitrarily on an equator (Fig. 3). 
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When (M, 7) = 0, the Hamiltonian and the additional integral are 




+ ^(73) 



(11) 



n 



F=Ml + {l- 1 l)Y J 



There is an arbitrary function t/(7 3 ) in This function means addition of an 

arbitrary "central" field with the center on a perpendicular to the plane of the 
Hookian potentials (Fig. 3). For example, one more Hookian center can be placed 
at the pole. This implies (see fl]) that the spatial problem of a point moving on a 
three dimensional sphere S 3 under the action of n Hookian centers on its equator is 
also integrable. 

Note that a Euclidean analogue of the problem in question is trivial, as even 
in Cartesian coordinates it yields n linear oscillators. In this case, the Hookian 
centers can be arbitrarily scattered in M. 2 . In the case of curved space, even on a 
two-dimensional sphere, the problem of motion in the field of three arbitrarily placed 
Hookian centers is not integrable, as simulations reveal chaos in this system. The 
quadratic integral F in (fTTj) is due to the fact that the problem is separable in the 
spherical coordinates (9, if). Indeed, the Hamiltonian H can be written as 



where 9, ip are the coordinates of the moving mass point, while <p i defines the position 




(12) 
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of the 2-th Hookian center on the equator (Fig. 3). The expression in square brackets 
is an additional integral of motion (fTTj) . 
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Figure 1: Location of real and "imaginary" centers on a plane 
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